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Q\ . Abstract: Atomic-like systems in which electronic motion is two dimensional 

are now realizable as "quantum dots". In place of the attraction of a nucleus 
there is a confining potential, usually assumed to be quadratic. Additionally, a 
perpendicular magnetic field B may be present. We review some recent rigorous 
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results for these systems. We have shown that a Thomas-Fermi type theory for 
the ground state is asymptotically correct when N and B tend to infinity. There 
are several mathematically and physically novel features. 1. The derivation of 
the appropriate Lieb-Thirring inequality requires some added effort. 2. When 
B is appropriately large the TF "kinetic energy" term disappears and a peculiar 



"classical" continuum electrostatic theory emerges. This is a two dimensional 
problem, but with a three dimensional Coulomb potential. 3. Corresponding to 
this continuum theory is a discrete "classical" electrostatic theory. The former 
provides an upper bound and the latter a lower bound to the true quantum energy; 
the problem of relating the two classical energies offers an amusing exercise in 
electrostatics. 
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I. INTRODUCTION 

Recent developments in semiconductor technology have produced remarkably small 
pieces of matter that behave like two-dimensional atoms. These "dots" are part of the 
world of mesoscopic physics which contains objects that are large compared to atoms but 
small enough to demonstrate typical quantum effects such as discrete energy levels and 
interference phenomena. 

The simplest quantum dot consists of an insulator, e.g. AlGaAs, and a semiconductor, 
e.g. GaAs. Between these two materials is a potential difference that confines injected 
electrons to a thin layer at the interface. The thickness of the layer ranges from 30 to 100 
angstroms while the diameter of the dot, typically of the order 300 to 8000 angstroms, 
is controlled by a metal gate on top of the insulator, or by edging a pattern through the 
insulator and the interface. (Recall that aBohr : = ^ 2 / (me 2 ), the radius of a hydrogen atom, 
is about half an angstrom, i.e. 0.5 x 10~ 8 cm). The number of injected electrons, N can 
conveniently be from zero to several thousand. 

Perpendicular to the dot a uniform magnetic field, £?, can be introduced, thereby 
making it possible to investigate the quantum mechanics of motion in a magnetic field. 
The quantized Landau levels can then become prominent. 

Although this system is very quantum mechanical, its parameters are sufficiently 
different from ordinary atoms and molecules that it is possible to demonstrate in the 
laboratory effects that would otherwise require the magnetic field of a neutron star (10 9 — 
10 12 gauss). 

• Energies are typically in units of meV (millielectron volts), not eV. 

• Large magnetic fields need only be several teslas (one tesla = 10 4 gauss). By a "large" 
field is meant one for which the effective Landau radius is comparable to the effective 
Bohr radius or smaller. 
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• The system is effectively two-dimensional because the electronic motion in the direc- 
tion perpendicular to the sandwich is essentially "frozen out" . The thickness, Lj_, is so 
small that each electron remains in the ground state, <j>(x±) = (\ / 2/L±) sin(7r;r_i_/L_i_). 
It costs a large kinetic energy to promote an electron to the next state 
(V2/L±) sin(27rx ± /L ± ). 

Whether or not such transitions can be totally disregarded remains to be proved, 
but this is the customary assumption and we shall make it here. From the mathematical 
point of view, the assumption of two-dimensional motion is the most interesting feature 
of this model. Some of our results here contrast markedly with those in our earlier work 
on conventional atoms in large magnetic fields [LSY1-3]. In this paper we present some 
recent work; details will appear elsewhere. 

Our iV-body system can be modeled by the following Hamiltonian, H^. 

N 

H N = Y J H] + el Yl (1-1) 

3 = 1 l<i<j<N 

with Xi G R 2 and where H 1 is the one-body Hamiltonian 

H 1 = -— — fiV -±a)\ 9 J^)s*B-(^)(i-^)b + V(x), (1.2) 
2m* V he J y \2m*cJ \2m*c J \ 2 J v h v ' 

with the following notation. 

The magnetic vector potential is A = ^B(—y,x). 

The potential V(x) is confining, which is to say that V(x) — > oo as |x| — > oo, and it 
is assumed to be continuous. For convenience in computations it is usually taken to be 
of the form V(x) = ^m*u 2 x 2 , but this is somewhat phenomenological and arbitrary. It 
is created by an external gate and it is therefore adjustable to a certain extent. This V 
replaces the Coulomb attraction to the nucleus, —Z/\x\, in a conventional atom. We write 
the potential generally as 

V(x) = Kv(x) 
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and let the coupling constant K and the magnetic field be N dependent, 



K = kN and B = (3N 

with (3 and k fixed. In this way the ground state energy, E 1 ^, will scale as N 2 as N — > oo. 
(This contrasts with iV 7 / 3 for a natural atom.) 

The operator S z is the perpendicular component of the spin and takes the values ±1/2 
for each electron. It commutes with H. 

The constant term in (1.2), — (fte/(2m*c))(l — \g*\/2)B, is included in order that the 
"kinetic energy" operator, H 1 — V(x), has a spectrum starting at zero. 

The Hilbert space is that appropriate for fermions, the antisymmetric tensor product 

JV 

A£ 2 (R 2 ;C 2 ). 
i 

The numbers m#, e*, g* are "effective" values that differ from the usual ones because 
of the electron's interaction with the crystal in which it finds itself. The values for GaAs 
are 

m* = effective mass = 0.07 m e i ec tron 

e* = effective charge = e- (dielectric constant) -1 / 2 = 0.3 e 
= -.03 

Corresponding to these values there are respectively a natural length, energy and magnetic 
field strength: 

a* = ft 2 /m*e 2 = 185 aBohr = 10~ 6 cm 
= ejm*/h 2 = 9 x 10" 4 Ry = 12 meV 

B m = eimlc/eh 3 = 6.7 x 10 4 G = 6.7 T 
Henceforth, we choose units such that fi = e* = 2m* = 1. The energy unit is then and 
if we moreover measure the magnetic field in units of the numerical factors in (1.1) 
and (1.2), except g*, can be dropped. The eigenvalues of (iV — A) 2 are the Landau levels 
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B, 3B, 5B in these units. The g-factor for normal electrons is 2, which is what appears in 
the Pauli operator, and which leads to eigenvalues 0, 2B, 4B, 6B, ... of (iV — A) 2 + 2S Z B, 
with the first being singly degenerate and the others doubly degenerate (as far as spin is 
concerned). In our case g* ^ 2. In fact, it is slightly negative in GaAs! This merely makes 
the bookkeeping a bit more complicated, but has no essential effect. A derivation of the 
effective g-factor in a solid is given in [KC, pp. 282-286]. 

We draw attention to the fact that the Coulomb repulsion in (1.1) is the three- 
dimensional potential, even though our system is two-dimensional. The absence 
of Newton's theorem makes our job a little harder. 

If the repulsion is omitted, and if V is quadratic, V(x) = \uj 2 x 2 : the problem is 
solvable. The spectrum of H 1 was determined by Fock [FV] in 1928, two years before 
Landau's paper. For (zV — A) 2 + \uj 2 x 2 the spectrum is given by 

E = ( ni - n 2 )B + (m +n 2 + l)[u 2 + B 2 ] 1 ' 2 (1.3) 

with 77.1,77-2 = 0, 1,2,.. .. It is remarkable that the simple spectrum, (1.3), gives a qualita- 
tively good fit to some of the data [KM, KLS]. 

One of the most important physically measurable quantities is the spectrum of H. 
Here we concentrate on the ground state energy and will demonstrate that the Thomas- 
Fermi (TF) energy, which scales precisely as iV 2 , is asymptotically exact as N — > 00 with 
k and (3 fixed. The iV — > 00 limit of turns out to be uniform in ft, i.e. given any e > 
we can find iV large enough, but independent of /?, so that A^~ 2 |i? TF — E®\ < e. 
II. THOMAS-FERMI THEORY 

TF theory is defined by the following functional, £ TF , of a nonnegative "density" 
function p = R 2 — > [0, 00). 
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£ TF (P) = Jj B (p(x))dx + J V(x)p(x)dx + D(p, p) 



R 2 



R 2 



with 



D{f,9) = \J f f(x)g(y)\x-y\- 1 dxdy. 



(2.1) 



(2.2) 



R 2 R 2 



The function j B (p) imitates the kinetic energy, (zV — A) 2 + g*S z B — (1 — \g*\/2)B. It is the 
smallest kinetic energy density for electrons in a 2-dimensional box with particle density p. 
For simplicity we give explicit formulas for spinless fermions. For electrons (with spin) this 
would correspond to having g* = (which is close to the value for GaAs) except for the 
trivial fact that one must replace j B (p) by 2j B (p/2) due to the spin degeneracy. Everything 
goes through, however, with obvious modifications, for electrons with g* ^ 0. If spin is 
ignored we have j B (0) = and j B (p) has slope 2nB for n(B/2n) < p < (n + l)(£?/27r), 
with n = 0, 1, 2, — See the graph. 



6B 

Ji 



3B 

Ji 



Kinetic energy density 
as a function of particle 
density, p . 




B/2JI B/Jt 3B/27C 2B/7C 
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W(x - y) := [\x - y\ 2 + 5Z 2 } 1 / 2 - [\x - y\ 2 + 4d 2 ] 1/2 , (2.3) 

where Sz ~ 100 angstroms is the thickness of the dot and d ~ 1000 angstroms is the 
distance to an external gate. This W has a short-distance and a long-distance cutoff. It 
is, however, positive definite as a kernel. Our methods can easily be adapted to prove the 
analogues of Theorems 1 and 2 with \x — y\~ x replaced by W everywhere. 
The TF energy is defined with V(x) = Kv(x), K = kN and B = f3N by 



E tf (k, p, N) = inf |£ TF (P) : p(x) > 0, J p = N, K = kN, B = (3N j (2.3) 
and it is evident that 

E TF ( K ,(3,N) = N 2 e TF (^(3) and p TF (x) = N pJ F (x) (2.4) 

where P\ and e TF are p and E TF in the special case N = 1. 

The functional £ TF is convex in p and it is strictly convex since the kernel \x — y\~ x 
is strictly convex. Standard methods yield 

1. THEOREM. For each N > and any g* there is a unique minimizer, i.e., a p 
such that S TF (p) = N 2 e TF {n, (5) and j p = N. 

The only novel feature here is the existence of a minimizer for all N, which is a 
consequence of the fact that V(x) — > oo as \x\ — > oo instead of V(x) — > as \x\ — > oo. 

The minimizer satisfies a TF equation, which we illustrate for the particular js defined 
above. Since j B (p) is not differentiable at the points p = n(B /2ir), we first have to define 
the set valued function for p > and n > given by 

j' B (p) = {2nB} if n(B/2n) <p<(n+ 1)(5/2tt) 

(2.5) 

j' B (p) = [2nB,(2n + 2)/B] if p = (n + 1)(B /2n) 
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where {x} denotes a single point x and [a, b] is the closed interval. The TF equation can 
now be written as 

-*,(,)-w-'.p + m{!&<**» «£>>J. p..) 

This holds also for arbitrary with a suitably modified j' B . The quantity /U is the physical 
chemical potential, i.e. 

» - a - E r 1 ("> 

with K and £? fixed. 

Not only is it true that p satisfies (2.6), but it is also true that any (nonnegative) 
solution pair (/z, p) to (2.6) is a minimizing pair, i.e. p is necessarily a minimizer for £ TF 
for some N and p is the chemical potential (2.7) at that N. Previously, such an assertion 
was proved [LS] for j a continuously differentiable function. The proof when j' is not C 1 
is a bit trickier and was carried out by Lieb and Loss [LL]. 

Our main theorem is the following 

2. THEOREM (Convergence of energy and density). Fix k, (3 andv(x). Then 
the quantum and TF energies converge in the sense that 

N- 2 \E tf (k,[3,N) - E q (k,/3,N)\ -> 

as N — > oo. The convergence is uniform in (3. The densities also converge as N — > oo: 

N~ 1 p Q (x) ^ pJ F (x) 

in the weak L 1 sense. 

Customary techniques, such as coherent states based on Landau orbitals [LSY2,3], 
are employed to prove this theorem, but they are insufficient for the following important 
reason. The function j B in (2.1) is zero in the interval < p < B/2n. Correspondingly, 



the quantum-mechanical kinetic energy operator (iV — A) 2 + g*S z B — (1 — \gJ/2)B has a 
whole band of zero modes. (Note: It makes no difference if g* ^ because that will merely 
entail a shift of the chemical potential.) These zero modes cause several problems which 
we turn to in Sect. III. 

Before doing so, however, let us comment on several unusual features of the TF density 
p(x). For one thing, it always has compact support. This is perhaps not surprising since 
the usual atomic TF density has compact support except when N = Z. A more interesting 
feature is that when B is large, < p(x) < B/2n for all x and hence j B (p(x)) = 0. This is 
not surprising. When B decreases, however, p(x) develops a flat spot because (2.6) causes 
p(x) to be anchored at B/2n for a set of positive measure in R 2 . The first situation for 
v(x) = x 2 , k = 10~ 4 and B = 8 tesla (with N = 50, so (3 = 6 x 10~ 3 in our units) is 
shown in Fig. 1. The upper part of the picture shows the density and the lower part the 
effective potential V e fi = V + * p (lower part) are shown. The second situation (B = 
7 tesla) is shown in Fig. 2. As B is decreased further the second Landau level starts to be 
occupied and the flat spot, caused by the anchoring at B, moves further out (Fig. 3, B = 5 
tesla). For a still lower field two flat spots appear because 2np is anchored at B and at 2B. 
Eventually more and more flat spots appear. Fig. 4 shows three flat spots at 2 tesla. From 
the point of view of the Schrodinger kinetic energy, we are seeing the filling of successive 
Landau bands. Also shown on these figures, for comparison, is a dotted curve representing 
the solution to a Hartree calculation of the quantum density for = 50. The curves for 
the TF theory were computed by Kristinn Johnsen while the Hartree calculations are due 
to Vidar Gudmundsson (see [PGM] for details about Hartree calculations for quantum 
dots). 
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III. TWO-DIMENSIONAL LIEB-THIRRING INEQUALITY WITH A 
MAGNETIC FIELD 

Here we consider spinless fermions to simplify the discussion. Given any normalized 
fermionic wave function ip(xi, . . . , xjv), with Xi E H d we can form the density 

Pip(x) = N J \ip(x ll x 2 , ■ ■ ■ ,x N )\ 2 dx 2 ■ ■ -dx N , (3.1) 

Rd(iV-l) 

so that / p^(x)dx = N. The Lieb-Thirring inequality relates T^, the kinetic energy of ip 
defined by 



N \ 

^V-Af-B^\ (3.2) 
i=i I 



to some integral of p^. If B = the original bound was > Kd J P^ +2 ^ d \ with Ka being 
a universal constant that is positive for all d > 1. 

If we omit the — B term in (3.2) then this original inequality continues to hold for 
all vector potentials A. If the — B term is retained, however, the situation is much more 
delicate because of the possibility of zero modes (which certainly exist for d = 2, for 
example, when B{x) has nonzero total flux). For d = 3 we proved [LSY3] that for some 
Ci > 

T^>d J jf{p^x))dx (3.3) 

where jg\p) is the function that appears in the d = 3 TF theory with a constant magnetic 
field. 

In two dimensions we would like to have a similar inequality, but the fact that j B (p) = 
for p < B/2tt in two-dimensions causes difficulties. The following is adequate, however. 

3. LEMMA (d = 2 kinetic energy inequality). WithT^ and p^ defined as above, 
with B constant in space and with js the kinetic energy function defined in Sect. II (i.e., 
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the energy per unit volume of spinless particles at density p), there is a constant C e > 0, 
for every < e < 1, such that 



Inequality (3.4) is equivalent to an assertion about the negative eigenvalues e\ < e2 < 
. . . < of H B = (zV — A) 2 — B — U(x) for an arbitrary nonnegative potential U{x). One 
shows that 



from which it follows that > f w(p^(x))dx with w being the Legendre transform of /. 
We are able to show that for every < A < 1 



and (3.4) follows easily from this. 

Our proof uses, as usual, the Birman-Schwinger kernel t/" 1 / 2 [(zV — A) 2 — B + E^U 1 ^ 2 . 
The hard part is to control the zero modes oih = (iV — A) 2 — B, i.e., the lowest Landau 
level. Two-dimensions is very different from three where the energy in the ^-direction 
precludes having a macroscopically degenerate level of zero energy states. We are able to 
isolate this lowest level with the aid of Ky Fan's inequality on singular values of sums of 
operators; the effect of this level is to produce the B J U term in (3.5). 

Technicalities aside, (3.4) presents a nontrivial mathematical problem. can, indeed, 
be zero; this is easily accomplished, even for positive density, by making use of all the 
lowest Landau level eigenfunctions. If we want to bound from below by the integral 
of a nonnegative, but nontrivial function of p, such as j B (p), we have to be sure that the 
function is zero for low density. The function j B , which is zero for < p < B/2ix, precisely 




(3.4) 




.x 




(3.5) 
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accomodates the maximum density in the lowest Landau band that is allowed by the Pauli 
principle. 

IV. EXCHANGE-CORRELATION ENERGY AND THE TWO CLASSICAL 
ELECTROSTATICS PROBLEMS 

The exchange-correlation energy of a normalized state if) is defined to be 



E,A4>) = ( 4> 



l<i<j<N 



^)-D(p^,p^), (4.1) 



with p^ defined as in (3.1) (but with a sum over spin states) and D{f,g) defined as in 
(3.2). By imitating the original proof in [LEI], i.e., using the "maximal function", we can 
prove 

4. LEMMA (Exchange inequality in 2 dimensions). There is a constant C > 
such that 

E ex (ijj) >-C f P^xf^alx. (4.2) 

This inequality is quite adequate for small and moderate B = (3N, for the simple 
reason that it can be controlled by / Jb{P^) as follows. 



P^<B/7T Ptp>B/tt 

1/2 



<Ny/Bfr+\ / P%\ {/p^} <|yf+ (const.) ^1 AT 3 / 2 . (4.3) 

The last line is 0(N 3 / 2 ) in the ground state, which is small compared to the total energy 
(which is 0(N 2 )) provided one can prove the "obvious" fact that = 0(N 2 ); this is easy 
to do. 
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For large f3 (4.3) is not optimal because it asserts that the error is /^^iV 3 / 2 . While 
this is 0(iV -1 / 2 ) relative to E®, it is not uniform in f3 and it is certainly not optimal. 
The source of the difficulty, once again, is the zero modes of the kinetic energy operator. 
We can do better by introducing the Classical Particle Energy, E^, defined as follows. 
First define the function on R 27V by 

N 

£^{xu...,x N ) = Ky^ j v(x i )+ ^ \xi-Xj\~ 1 . (4.4) 

i=l l<i<j<N 

Then E^ is defined by 

E*{N) = inf{££(xi, . . . , x N ) : x % G R 2 }. (4.5) 
Corresponding to this particle energy is the Classical Continuous Energy: 

£%{p) = K J p(x)v(x)dx + D(p, P) (4.6) 

R 2 

E°{N) = mf{£%(p) : p(x) >0,Jp = N}. (4.7) 

The classical continuous problem is the same as the TF problem, but without the 
kinetic energy j B (p) term and hence, in particular, it does not involve B. As a result of 
the lack of a f p 2 control, the minimizer, which one can show always exists (since V(x) — > oo 
as \x\ — > oo) and is unique, might turn out to be a measure! Since V is continuous, this 
makes sense, but we question whether this actually happens. (Note: There are Hausdorff 
measures in two-dimensions that have finite energy D{p,p).) However, as discussed just 
before Lemma 7, measures to not occur when V G C ,a . 

The classical continuous problem has the same simple scaling property as the TF 
problem, namely 

E%(N) = N 2 E^(1) and p N (x) = N Pl (x), (4.8) 
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where Pat is the 



(with k fixed) for J p = N and K = kN. 



The minimization problem (4.7) leads, of course, to the following Euler-Lagrange 
equation for p that parallels the TF equation (2.6). For N = 1, K = k 



with p being the chemical potential. [Note: if p is a measure the left side of (4.9) should 
be > p everywhere and = p p-almost everywhere. 

Eq. (4.9) states that the potential generated by p plus the confining potential V add 
up to a constant when p > 0. If we were in three-dimensions instead of two, the solution 
is easy. Just apply —A to both sides. If v(x) = x 2 one has that p is a constant, C, inside 
some ball of radius R (with C = 3k/(2tt), 2kR 3 = 1, p = 1/R + kR 2 ). 

The solution to (4.9) in two-dimensions is much more difficult. However, when v(x) = 
x 2 it can be solved: 



for suitable constants C and R. The computation of the convolution (right side of 

(4.10)) is an amusing exercise! 

We have three non-quantum energies before us, E tf (k, ft, N), E^-(N), E^-(N). In 
order to discuss their relationship we first need a technical lemma 

5. LEMMA (Finite radius). Fix k and set K = kN. 

There is a radius R v (k), independent of N, such that for every N the infimum for 
can be sought among densities (in the case of E TF and E c ), or particle 
coordinates (in the case of E p ) that have support in the ball of radius R v (k) centered at 
the origin in R 2 . 

The proof of Lemma 5 is fairly simple. If there is some charge outside some large ball, 
R v {k), then there is some location inside the ball of unit radius to which this charge can 




(4.9) 




(4.10) 
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be moved in such a way that the energy is lowered. Here, we make use of the fact that 
V(x) — > oo as \x\ — > oo. 

In the remainder of this section we discuss the relationship between the two classical 
energies E^-(N) and E^-(N). In the next section we discuss the relationship of these 
classical energies to E TF . We first point out that Lemma 4 cannot be used in order to 
compare E^ and E%. In fact, what corresponds to the minimizing "density" p^ for the 
problem is a sum of delta functions and thus f p 3 ^ 2 is meaningless (if we approximate 
the delta functions by continuous functions this error term will go to infinity). 

Using an electrostatics inequality of Lieb and Yau [LY] we are able to prove that 
Ej^-(N) is indeed a good approximation to the particle energy Ej^(N). That electrostatics 
inequality generalizes (4.2) in a certain sense. It is the following 

6. LEMMA (The interaction of points and densities). Given points x±, . . . , xn 
in R 2 , we define Voronoi cells T\, . . . , Tjy C R 2 by 

Tj = {y E R 2 : \y - xj\ <\y-x k \ for all k ^ j} . 

These Tj have disjoint interiors and their union covers R 2 . We also define Rj to be the 
distance from Xj to the boundary ofTj, i.e., Rj is half the distance of Xj to its nearest 
neighbor. Let p be any (not necessarily positive) function on R 2 . (In general, p can be 
replaced by a measure, but it is not necessary for us to do so.) Then 

N f 

\ x i~ x j\ 1 >-d(p,p) + J2 p(y)\y~ x j\ ld v 

l<i<j<N j = l ^ R2 

N N . 

+ 1 sJ2 R J 1 -J2 Piy^y-x^dy . (4.11) 

3=1 3=1 Tj 

This comparison of E^(N) and E^-(N) requires that the minimizing density for the 
continuous energy has some regularity. A-priori the minimizing density for could 
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be a singular measure and if this is the case we cannot expect to get a good estimate of the 
difference between and E^. To ensure that the minimizer for is neither a measure 
not too singular a function we need the confining potential V to be sufficiently regular. 
Although we do not believe this to be optimal, we find it convenient to assume that V is 
in C 1,a (Holder continuous first derivative with exponent a > 0). We prove that p is then 
in L 9 (R 2 ) for some exponent q > 2 that depends on a. With this in hand, we are then 
able to prove 

7. LEMMA (Classical correlation estimate). Assume that V is a potential in 
C 1,a . For all N we have 

E%(N) - b(K)N 3 / 2 < E%(N) < E%(N) - a(^)A 3 / 2 , (4.12) 

where cl(k) and &(«) are positive constants independent of N that depend on the radius 
R v (k) of Lemma 5 and on the exponent q (and its dual p = q/ (q — 1) ). 

a{n) = (8RM)- 1 

with Pi being the minimizer for £^ for N = 1, i.e., the minimizer for . 

V. OUTLINE OF THE PROOF OF THE LIMIT THEOREM 

We now give a very sketchy outline of the proof of Theorem 2. We shall mainly 
emphasize the steps in the proof that do not use standard techniques. In particular, we 
shall stress the use of the estimates given in the two previous sections. 

The aim is to prove that for a given e > we can find an electron number N(e, k), 
depending on e and k but independent of (3, such that > N(s, k) implies 

N- 2 \E tf (k,(3,N)-E®(k,(3,N)\ <e (5.1) 
16 



for all (3. To prove this we have to treat small and large (3 differently. We therefore choose 
P £ (to be specified below) and consider (3 < (3 £ and (3 > (3 e . In both cases the proof 
consists of finding upper and lower bounds to E®(k, (3, A). As usual, the upper bound is 
a relatively simple variational argument. In both cases we construct a trial wave function 
using magnetic coherent states similar to the ones used in [YJ] and [LSY1-3]. The lower 
bounds are more difficult; the case of large (3, in particular, requires new ideas — again due 
to the presence of a large number of zero modes. We discuss this large (3 case first. 

Case 1. (Lower bound on E® for (3 > (3 e ): By simply ignoring the kinetic en- 
ergy operator, which we had normalized to be positive, we have the obvious inequality 
EQ(k,P,N)>E£(N). 

From Lemma 7 we can therefore conclude that 

E Q (n,f3,N) > E%(N) -b(K)N 3/2 . (5.2) 

The final step in the proof is to relate E^-(N) to E tf (k, f3, N). In order to do this 
we choose a bounded density Pi with J p 1 = 1 such that £%{Pi) < E^(l) +e/2. Note that 
does not depend on f3. We may then assume that (3 e is so large that P\{x) < /3 £ /2n 
for all x. For all (3 > (3 e we therefore have jp^{Np) = (for simplicity we are here again 
ignoring spin). 

Consequently, using p = NPi as a trial density in the TF functional we find that 

£ TF (^,/3, A) < £ TF (p) = £g(p) = N 2 [E c K {l)+e/2] < E°{N) + |a 2 . (5.3) 

Comparing (5.2) and (5.3) we immediately see that N(e, k) can be chosen so that (5.1) 
holds. 

In the special case of v (x) = x 2 we know that the minimizer P\ of 8% given by (4.10) 
is continuous and bounded. In this case we can therefore choose P\ = P\ and we can make 
the explicit choice (3 e = sup x 2nPi(x) = 2nPi(0) = CR. 
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Case 2. (Lower bound on E Q for f3 < f3 e ): In this case we follow the standard route of 
using inequalities (4.2) and (4.3) to reduce the many-body problem to a one-body problem. 
Combining (1.1), (4.1-3) we obtain 

JV 



\H N m = + D(p+, p+) - C(k, e)N 3 / 2 , 



(5.4) 



where C(k, e) is chosen as an estimate on the factor in { } on the right side of (4.3), i.e. 



y/Pe/ir + (const.) ^/T^/N < C(k,e). 

To relate (5.4) to the Thomas- Fermi problem we use the inequality 



< D{p^ - p TF , - p rt ) = - ( V 



JV 



3 = 1 



,TF 



* X, 



^)+D(^,^)+D(p TF ,p TF ) 



which is a consequence of the positive definiteness of the kernel \x — y\~ x . Inserting this 
in (5.4) gives 

JV / JV 



\H 



N\ 



>E ^l+E^^N" 1 



$ ) -D(p TF ,p TF ) -C( K ,e)N s '' z . (5.5) 



3/2 



The importance of (5.5) is that the quantum energy, (iJ)\Hn\iP), is bounded below by the 
expectation of the one-body operator H 1 + p TF * 

To prove Theorem 2 it remains to do a semiclassical analysis (with iV -1 / 2 playing 
the role of an effective Planck's constant) of the sum of the negative eigenvalues of this 
one-body operator. This semiclassical analysis again uses the magnetic coherent states 
and is very similar to the analysis in [LSY3]. 



Acknowledgements We thank Kristinn Johnsen and Dr. Vidar Gudmundsson for valu- 
able comments and for allowing us to use the pictures of densities and effective potentials 
that they prepared. 



18 



REFERENCES 

[FV] V. Fock, Bemerkung zur Quantelung des harmonischen Oszillators in Magnetfeld, 

Z. Phys. 47, 446-448 (1928). 
[KC] C. Kittel, Quantum Theory of Solids, Wiley (1963). 
[KLS] A. Kumar, S.E. Laux and F. Stern, Electron states in a GaAs quantum dot in a 

magnetic field, Phys. Rev. B. 42, 5166-5175 (1990). 
[KM] M.A. Kastner, Artificial atoms, Physics Today 46, 24-31 (1993). 
[LEI] E.H. Lieb, A lower bound for Coulomb energies, Phys. Lett. 70A, 444-446 (1979). 
[LL] E.H. Lieb and M. Loss, unpublished section of a book on stability of matter. 
[LS] E.H. Lieb and B. Simon, Thomas-Fermi theory of atoms, molecules and solids, Adv. 
in Math. 23, 22-116 (1977). 
[LSY1] E.H. Lieb, J. P. Solovej and J. Yngvason, Heavy atoms in the strong magnetic field of 

a neutron star, Phys. Rev. Lett. 69, 749-752 (1992). 
[LSY2] E.H. Lieb, J. P. Solovej and J. Yngvason, Asymptotics of heavy atoms in high magnetic 

fields: I. Lowest Landau band region, Commun. Pure Appl. Math, (in press). 
[LSY3] E.H. Lieb, J. P. Solovej and J. Yngvason, Asymptotics of heavy atoms in high magnetic 
fields: II. Semiclassical regions, Commun. Math. Phys. 161, 77-124 (1994). 
[LY] E.H. Lieb and H.-T. Yau, The stability and instability of relativistic matter, Commun. 

Math. Phys. 118, 177-213 (1988). 
[M] P.L. McEuen, E.B. Foxman, J. Kinaret, U. Meirav, M.A. Kastner, N.S. Wingreen and 
S.J. Wind, Self consistent addition spectrum of a Coulomb, island in the quantum Hall 
regime, Phys. Rev. B, 45, 11419-11422 (1992). Cf. eqn. (3). 
[PGM] D. Pfannkuche, V. Gudmundsson, P.A. Maksym, Comparison of a Hartree, a Hartree- 
Fock, and an exact treatment of quantum dot helium, Phys. Rev. B, 47, 2244-2250 
(1993). 

[YJ] J. Yngvason, Thomas-Fermi theory for matter in a magnetic field as a limit of quantum 
mechanics, Lett. Math. Phys. 22, 107-117 (1991). 



19 




20 



